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ABSTRACT 

We present the first numerical simulations that self-consistently follow the formation 
of dense molecular clouds in colliding flows. Our calculations include a time-dependent 
model for the H2 and CO chemistry that runs alongside a detailed treatment of the 
dominant heating and cooling processes in the ISM. We adopt initial conditions char- 
acteristic of the warm neutral medium and study two different flow velocities - a slow 
flow with 'L'flow = 6.8 km s~^ and a fast flow with v^^^ = 13.6 km s~^. The clouds 
formed by the collision of these flows form stars, with star formation beginning after 
16 Myr in the case of the slower flow, but after only 4.4 Myr in the case of the faster 
flow. In both flows, the formation of CO-dominated regions occurs only around 2 Myr 
before the onset of star formation. Prior to this, the clouds produce very little emission 
in the J = 1 ^ transition line of CO, and would probably not be identified as molec- 
ular clouds in observational surveys. In contrast, our models show that H2-dominated 
regions can form much earlier, with the timing depending on the details of the flow. In 
the case of the slow flow, small pockets of gas become fully molecular around 10 Myr 
before star formation begins, while in the fast flow, the first Il2-dominated regions oc- 
cur around 3 Myr before the first prestellar cores form. Our results are consistent with 
models of molecular cloud formation in which the clouds are dominated by "dark" 
molecular gas for a considerable proportion of their assembly history. 
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1 INTRODUCTION 

The life-cycle of molecular clouds sets the time-frame over 
which the star formation process can occur. As such, molec- 
ular cloud formation has received much attention in recent 
years, as we try to distinguish between theories for 'rapid' or 
'slow' star formation. In a simplified picture, the cycle can 
be thought of as comprising three distinct phases: assembly 
of the cloud, the formation of stars within the cloud, and the 
cloud's eventual dispersal (often linked to the termination 
of the star formation process). While we have made some 
progress on understanding the latter two stages of the cloud 
life-cycle (see e.g. |Tamburro et al.|2008 Jeffries et al.|20lT ), 
the first stage is more uncertain, and has prompted much 
debate (see e.g. Ballesteros-Paredes, Hartmann & Vazquez- 



Semadeni"1999 ; Elmegreen 2000; Hartmann et al. 2001 ; Tas- 
sis Mouschovias 2004; Mouschovias, Tassis Kunz 2006 
Blitz et al.|[2007| |Elmegreen||2007| [Tamburro etUT 



Pagani, Roueff Lesaffre||2Qll ) 



2008^ 



The idea of the colliding flow model for the formation 
of giant molecular clouds (GMCs) offers a simplified pic- 
ture of how clouds can be assembled: two warm, marginally 
supersonic flows of gas collide head-on, and a dense molec- 



ular cloud builds up in the shocked layer ("Elmegreen"" 1993 



Walder FQliiiIl|1998a|bl [Hennebelle Perault.1999, .2000 



&: Hen nebelle'' 2005f Vazquez-Sem adeni et al.||2006| [Heitsch 
et al. 2006; H ennebelle et al.||200 8; Ban erjee et al.||2009| ). 



Koyama Inutsuka||2000[ |2002| [Bergin et al.||2004| [Audit 



This model has several attractive features. First, there is no 
need to invoke some ad hoc model for driving turbulence in 
the cloud, as the collision between the flows and the ther- 
mal instability that occurs within the shocked gas naturally 
yield a linewidth-size relationship comparable to those ob- 
served within Galac tic GMCs ([Koyama Sz Inutsuka 2002[ 
Bonnell et aI][20Q6l [Pobbs k Bonnell[[2QQ7| [Hennebelle etj 



al. 



2007 Heitsch, Hartmann & Burkert 2008). Second, to 



some zeroth approximation, the model probably captures 
how most molecular clouds form; or at least those in spi- 
ral galaxies, where the gas streams have ample opportu- 
nity to collide (Dobbs 2008). Furthermore, it seems that 



the timescale involved in the buildup of dense regions is 
reasonably insensitive to the strength of the magnetic field 
(Banerjee et al. 2QQ9J . Finally, the timescale for cloud for- 
mation, and the onset of star formation can then be iden- 
tified with the timescale for the onset of the gravitational 
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instability in the shocked layer. The model is very much in 
the spirit of 'fast' or 'dynamical' star formation, where the 
whole process of cloud assembly and gravitational collapse 
to form stars occurs on the dynamical timescale - or 'cross- 
ing time' - of the system ( [Ballesteros-Paredes, Hartmann 
Vazquez- Semadeni 1999; Elmegreen 2000; Hartman n et al. 
2001| [Vazqu ez- Semad eni et al.^^2007J , yielding groups of 
stars that have age-spreads smaller than their crossing times 



(see also the discussions in Mac Low & Klessen 2004 and 
McKee Ostriker|2007] ) . 

A crucial feature of the colliding flow model is the 
better-than-isothermal cooling that is offered by the inter- 
stellar medium as the gas is compressed. In particular, the 
cooling provided by via its 158 /xm fine-structure transi- 
tion permits the gas to undergo an isobaric thermal instabil- 
ity ( |Field||1965| ) that can produce a 'two-phase' interstellar 
medium (ISM), with warm diffuse gas in pressure equilib- 
rium with cold dense gas. In the conditions appropriate for 
the local ISM, the warm phase, known as the warm neutral 
medium ( WNM) , has a temperature of around 5000 K and a 
number density of around 1 cm~^, while the cold phase (the 
cold neutral medium, or CNM) has a temperature of around 
50 K and a number density of around 100 cm~^ (Wolfire et 
~ Ferriere|2001t . 



8tDU995^ 



2003 



An important issue that remains to be addressed is the 
ability of molecular gas to form in these flows. In order to 
form a GMC, it is necessary to convert a significant mass of 
hydrogen from atomic to molecular form, and to do this in 
a timescale that is comparable to or less than the assembly 
time of the GMC. However, with the exception of a few 



one-dimensional studies (e .g. [Hennebelle Perault 1999 



2000 Koyama &: Inutsuka 2000 Bergin et al. 2004), most 



previous simulations of colliding flows have not modelled 
the chemical evolution of the gas. Therefore, although they 
have demonstrated that large clouds of cold gas can be built 
up in this fashion, they have been unable to address what 
fraction of this gas will be molecular and what fraction will 
remain atomic. Also, they have have been unable to make 
solid observational predictions for what the cloud formation 
process should look like. 

The formation of H2 in the turbulent ISM has been 
modelled in a number of simulations that do not directly in- 
vestigate the colliding flow paradigm. For example, on large 



scales, the work of Dobbs et al. (2008| has demonstrated 
that 'molecular' clouds can be formed via compression of 
the gas in the spiral arms of galaxies. For scales greater 
than around 10 pc, this study demonstrated that H2 can 
form on timescales shorter than the crossing time of the 
spiral-arm passage. On smaller scales, j Glover Mac Low| 
( 2007b ) have shown that H2 can easily form on a dynamical 



timescale in gas with number densities of around 100 cm~^, 
thanks to the boost to the H2 formation rate provided by 
the many transient turbulent density compressions. Similar 
results were also recently reported by Micic et al. ( 2012| for 
the case of compressively-driven turbulent flows, and [Mac] 
Low & Glover (2012| have shown that the ratios of H2 to 
atomic hydrogen produced by these small-scale models af- 
ter a single gravitational free-fall time are in good agreement 
with the observed values (Blitz k Rosolowsky"2004"2006V If 
one combines the results of these studies, then it would seem 
clear that forming clouds of H2 is fairly easy, and that the 



colliding flow model should have no difficulty in producing 
large quantities of gravitationally-unstable molecular gas. 

However it is via CO, not H2, that molecular clouds are 
observationally defined, and so any formation model must 
necessarily produce enough CO to allow the cloud to be 
identified as a GMC by observational surveys. To date, the 
formation of CO has only been studied directly in mod- 



els of turbulent, pre- assembled clouds (see e.g. Glover et al. 
20T0I [Glover Mac Low||2011| [Glover Clark[[2012a|bT 



Shetty et al.[[2011a|b[ for models that treat both the hydro- 
dynamics and the chemistry, or [Rollig et al.f '^007 for an 
entry point into the extensive astrochemical literature deal- 
ing with static cloud models). The issues of how the CO 
formation timescale relates to the timescale on which the 
gas is assembled, or how it relates to the emergence of star- 
forming regions within the cloud have yet to be explored 
self-consistently. The closest that any previous colliding-fiow 
study has come to this is the work of [Heitsch k Hartmann[ 
( 2008[ hereafter HH08). They took a simplified approach to 
this problem, and suggested that any region in a colliding 
fiow that has a temperature of less than 50 K and a mean 
visual extinction of Ay = 1 or higher will have a high CO 
abundance. Using this approximation, they found that col- 
liding fiows could produce large, observable molecular clouds 
on timescales shorter than the star formation timescale (by 
around 4-5 Myr in the fiows that they looked at). However, 
given the non-equilibrium nature of the CO chemistry, and 
the complicated relationship between density, temperature 



and extinction found by Glover et al. (2010), it is unclear 



whether the assumptions made by HH08 are correct, and it 
is therefore worth revisiting this type of calculation with a 
more sophisticated treatment of the chemistry. 

With this in mind, we revisit the colliding-fiow model of 
cloud formation with a self-consistent time-dependent chem- 
ical network - coupled to a thermodynamical model of the 
ISM - that runs alongside the hydrodynamics. The chemi- 
cal network is able to evolve the abundances of key species 
such as H, H^, H2, C, C^ and CO, and allows us to follow 
the transition of the gas from the WNM to the CNM and 
onwards to the cold molecular cores that are the sites of star 
formation. 



2 SIMULATIONS 

2.1 The numerical model 

We model the evolution of the gas in our simulations using 
a modified version of the Gadget2 smoothed particle hy- 



drodynamics code (Springel 2005j). We have extended the 
publicly available code in several respects. First, we have 
added a sink particle implementation, based on the prescrip- 
tion in Bate, Bonnell & Price (1995), to allow us to follow 



the evolution of the gas beyond the point at which the first 
protostar forms. Our particular implementation is the one 
first described in Jappsen et al. (2005). In addition, we have 
implemented an external pressure term (e.g. Benz| 199*0 ) into 
the SPH equations that enables us to model a constant pres- 
sure boundary, as opposed to the vacuum, periodic, or 'fixed' 
boundary conditions that are the only choices available in 
the standard version of Gadget2. The details of this term 



can be found in Clark et al. (2011) 
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Figure 1. Column density images showing the evolution of the 6.8 km s"-*^ (slow) flow at 3 points before the formation of the flrst sink 
particle. The top images show the x-y plane, perpendicular to the flow, while the bottom panels show the cloud as seen looking down 
the aj-axis onto the z-y plane - that is, looking along the flow, onto the shock plane. 



Our treatment of the gas-phase chemistry combines the 
hydrogen chemistry network introduced in [Glover Mac| 
Low ( 2007a|b ) with the treatment of CO formation and de- 
struction proposed by Nelson & Langer (1999). This com- 



(Bohlin, Savage & Drake 1978 Draine & Bertoldi 1996), 



bined network was tested against other simplified chemical 
networks for CO formation and destruction by [Glover fc| 
Clark (2012b), who showed that it does a very good job 



of reproducing the results of more complex networks (e.g. 
Glover et al.||2010 ) for C^, C and CO, while incurring only 
one-third of the computational cost. Further details regard- 
ing our treatment of the chemistry can be found in [Glover fc| 
Clark (2012b). Note that the simulations presented in this 



paper do not include the effects of the freeze-out of CO onto 
dust grains. In Galactic star-forming clouds, CO freeze-out 
has been shown to have only a very minor effect on the 
thermal balance of the gas ( Goldsmith[[2001 ) . 



To model how the interstellar radiation field (ISRF) 
penetrates into the cloud, we make use of the TreeCol al- 



gorithm (Clark, Glover & Klessen 2012). This algorithm pro- 



vides each SPH particle with a Atv steradian column density 
map of the sky, stored in 48 Healpix pixels ("Gorski et al. 
[2005 ) . The column density information can then be used to 



calculate the attenuated ISRF at the location of each parti- 
cle - information that is vital for calculating the dust tem- 
peratures and photoelectric emission rates, as well as com- 
puting the photo-dissociation rates in the chemical networks. 
To convert from the column density to the dust extinction, 
we adopt the relationship Ay = 5.348 x 10~^^[A^H,tot/lcm~^] 



where A^H,tot is the column density of hydrogen nuclei. We 
also employ TreeCol to compute maps of the H2 and 
CO column densities, and use these to determine shield- 
ing factors due to H2 self-shielding, CO self-shielding and 
the shielding of CO by H2. For H2 self-shielding, we take 
the a shielding function given in Draine & Bertoldi (1996), 



while for the other two processes, we make use of the val- 
ues tabulated by Lee et al. (1996). We would not expect 



to obtain significantly different results were we to use the 
updated treatment of CO self-shielding presented in | Visser"^ 
van Dishoeck, k Black| (|2009|. 



2.2 Details of the colliding flow 

We start our flow at a number density of 1 cm~^ (note that 
we will refer to 'number density' simply as 'density' in what 
follows) and a temperature of 5000 K. The geometry of the 
flow is a rectangular cuboid of length zb 111.5 pc in the x- 
direction and a width of ± 36 pc in y and z. The total mass 
in this volume is 4 x 10^ M©. We apply a flow of +Vflow 
km s~^ for particles at negative x and —inflow for particles 
at positive x, such that the flows move along the x-axis 
towards zero. In this paper we look at two flow speeds, 6.8 
km s~^, and 13.6 km s~^, which we will refer to as the 'slow' 
and 'fast' flows respectively. The time taken for all of the 
gas to reach the centre is therefore ^16.6 Myr in the slow 
case, and ^8.3 Myr in the fast case. The isothermal sound 
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Figure 2. As FigurefT] but for the 13.6 km s ^ (fast) flow. Note that the output times examined here are different from those in FigurefT^ 



speed in the gas is initially 5.6 km s~^, so our flow velocities 
have initial Mach numbers of = 1.22 and A4 — 2.62, for 
the slow and fast flows, respectively. We also impose a low 
level of turbulence on the flows that has a power spectum of 
P(k) oc A:~^, with turbulent velocities scaled such that the 
rms turbulent velocity is given by Vrms = 0.2 Cg. 

We use 8,000,000 SPH particles to model the gas, and 
the particles are initially placed on a grid. When comput- 
ing fluid properties, we require that each particle has 70 
neighbours. The resulting mass resolution is therefore 0.7 
M© ( [Bate k Burkert|19 97). We have found in previous work 
that a mass resolution of around this value is suflficient 
for determining the star formation rate ( | Glover Clark 
2012a), but note that it is probably not sufficient to al 



low us to draw any strong conclusions regarding the stel- 
lar initial mass function (IMF). We distribute the particles 
on a uniform grid at the start of the simulation. We then 
evolve the flow, using a value of 1.5 for the Gadget2 vis- 
cosity parameter. We prevent immediate expansion of the 
flow through use of a confining external pressure, which we 
set to be Pext = 5000 /cb, where /cb is the Boltzmann con- 
stant. This value is comparable to the expected equilibrium 
pressure in the warm neutral gas given our choices for the 
metallicity, UV field strength and cosmic ray ionization rate 
described below. 

We consider solar metallicity gas, and adopt values for 
the initial ele mental abundances of c arbon, silicon and oxy- 
gen given by 



Sembach et al. 

^ and xo. 



(2000|: xc,tot = 1.4 x 10" 

4 



number of hydrogen nuclei. We assume that the carbon and 
silicon start in the form of and Si^, owing to photoion- 
ization by the ISRF, and that the oxygen starts in neutral 
atomic form. The initial abundance in the gas was set to 
Xh+ = 6.5 X 10~^, which is approximately the value yielded 
by equilibrium between cosmic ray ionization and radiative 
recombination. The remainder of the hydrogen was assumed 
to be in atomic form. The dust-to-gas ratio was taken to be 
0.01, the canonical value for the local ISM. For the dust 
opacities, we use two main sources of data. For the long 
wavelength opacities ( A > 1 /xm) , we use the values given by 
Ossenkopf & Henning ( 1994) for non-coagulated grains with 



thick ice mantles. For shorter wavelengths, we use the val- 



ues quoted in Mat his, Mezger Panagia] ( |l983 ) . In reality, 
dust grains in the warm ISM are highly unlikely to have ice 
mantles, but we note that in this regime, the particular val- 
ues chosen for the long wavelength grain opacities have very 
little influence on the thermal evolution of the gas. They be- 
come significant only in regions with high dust extinctions, 
where both gas and dust are typically very cold, and in these 
circumstances the grains will likely develop ice mantles. 

For the interstellar radiation field, we use the parame- 
terisation given by Black ( 1994) for optical and longer wave- 



lengths, and the Draine (1978) fit in the ultraviolet regime 



xsi,tot = 1.5 X 10 and a:o,tot = 3.2 x 10 , where Xi de- 
notes a fractional abundance, by number, relative to the 



The cosmic ray ionization rate of atomic hydrogen in all of 
our simulated clouds was chosen to be Ch = 3 x 10~^^ s~^. 
The cosmic ray ionisation rates for the other major chemical 
species tracked in our chemical model were assumed to have 
the same ratio relative to the rate for atomic hydrogen as 
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given in the UMIST99 chemical database ( |Le Teuff, Mihar 
k Markwick^,20QQJ . 

As we include the effects of self-gravity in our simula- 
tions, the gas is eventually able to form bound regions that 
start to undergo run-away gravitational collapse. In order 
to allow us to continue the simulations beyond the point 
at which the first of these regions collapses, we make use 
of sink particles. Regions that reach a density of 10^ cm~^ 
are considered as candidates for sink particle creation. To 
ensure that sink particles are not created out of transiently 
collapsing gas that should later re-expand, the SPH parti- 
cles at the centre of the region have to pass a series of tests 
to verify that they are actually gravitationally bound and 



collapsing, as outlined by Bate, Bonnell & Price (1995). In 



our study, the sink particles have an accretion radius of 600 
AU, and their gravitational interactions are softened on a 
scale of 100 AU. We also limit the distance between a new 
sink forming and a neighbouring sink to 1200 AU. Sink par- 
ticles are permitted to accrete SPH particles once the gas 
falls inside the accretion radius, provided that the candidate 
SPH particle is not only bound to the candidate sink par- 
ticle, but more bound to that sink than to any other. The 
SPH particle's mass and momentum are then transferred to 
the sink particle. In this manner, the mass in sink particles 
can be used as a proxy for the mass in stars. 

In addition to the two simulations discussed in the fol- 
lowing sections, we also performed two additional simula- 
tions that we do not discuss in detail here. The first of these 
was a lower resolution simulation with only 200,000 SPH 
particles that we used to check our numerical convergence. 
We found that while many of the features that we discuss in 
this paper were delayed with respect to the higher resolution 
runs, the basic physical and chemical traits were the same. 
The delay - which was around 1 million years - is expected, 
since when the resolution is poor, the small scale fluctua- 
tions are washed out, and the shock heating is spread over 
a larger mass of gas. The second additional simulation em- 
ployed a different realisation of the random turbulent field 
that we use to inject inhomogeneities into the flow. Again, 
we did not find any significant differences that would alter 
the conclusions that we draw below. 



3 EVOLUTION OF THE FLOWS 

In Figures [l] and [2] we show images of the column density 
of the gas, looking both perpendicular to the axis of the 
flow (upper panels) and along the axis (lower panels). Im- 
ages are shown for three output times: the time at which 
the first sink particle forms (hereafter referred to as tsF, de- 
noting the onset of star formation in the clouds), plus two 
earlier times that are approximately 50% and 75% of tsF- 
The images show a sequence of events similar to those seen 
in previous colliding flow models. In the region where the 
flows collide, a thermal instability is triggered which pro- 
duces a dense, cold layer of gas with significant turbulent 
sub-structure. Initially, the small-scale structure in this cold 
layer is not self-gravitating, but eventually some small dense 
regions become gravitationally bound and go into run-away 
gravitational collapse, leading to the formation of stars. In 
the case of the slow flow, star formation begins after roughly 
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Figure 3. Evolution with time of the maximum density (blue, 
solid line) and minimum temperature (red, dashed line) in the 
slow flow (top panel) and the fast flow (bottom panel). Note that 
at any given instant, the coldest SPH particle is not necessarily 
the densest, and so the lines plotted are strictly independent of 
one another. 



16 Myr, while the fast flow begins to form stars after only 
4.4 Myr. 

A more quantitative overview of the evolution of the 
flow is provided by an examination of the density and tem- 
perature evolution of the gas. As molecular clouds are typ- 
ically much denser and colder than our initial conditions, 
we first look at the evolution of the maximum density and 
minimum temperature in the gas, and how this develops as 
the flows come together. These quantities are plotted in Fig- 
ure |3] and were calculated for the densest and coldest SPH 
particles in the simulation at any point. Note that although 
dense gas tends to be cold, the densest particle is not nec- 
essarily the coldest particle, and so at any given time, the 
lines can be defined by different SPH particles. 

Initially, the two flows exhibit very similar behaviour: 
the minimum gas temperature drops to around 30 K within a 
few Myr, accompanied by a gradual rise in the peak density. 
The time taken for this initial phase of cooling is longer for 
the slow flow than for the fast flow. This can be understood 
as a consequence of the difference in the initial post-shock 
density in the two cases. The shocked gas in the collision re- 
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Figure 4. The plots on the left show the mass of gas that sits above a given density threshold as function of time, while those on the 
right show the fractions of the cloud mass that lie within given temperature ranges. The top panels are for the 6.8 km s"-*^ flow, and the 
bottom panels follow the 13.6 km s "-"^ flow. Note that the scales on the time axis differ between the top and the bottom panels. 



gion of the fast flow has a higher density, owing to the higher 
Mach number of the flow, and hence cools more rapidly than 
is the case with the slow flow. 

At first, the shocked gas has a higher thermal pressure 
than the unshocked gas surrounding it, although it is pre- 
vented from expanding in the x-direction by the ram pres- 
sure of the inflowing gas. As the gas cools, however, it first re- 
establishes thermal pressure equilibrium with its surround- 
ings and then maintains this equilibrium by increasing its 
density as its temperature decreases. By the time the gas has 
cooled to 30 K, its mean density has increased by roughly 
two orders of magnitude (see Figure|4|. The cold gas has also 
become turbulent, leading to a much larger increase in the 
peak density of the gas, owing to the influence of turbulent 
compressions occurring within the cold slab. 

The initial cooling phase is complete after around 4 Myr 
for the slow flow, and roughly 1.5 Myr for the fast flow. At 
this point, there is a marked change in the behaviour of the 
maximum density and minimum temperature. In the slow 
flow, we see that both quantities remain roughly constant 
for around 10 Myr. In the fast flow, a similar phase occurs 
for roughly 2.5 Myr, although in this case there is a slow 
but steady decrease in the minimum temperature and a sim- 
ilar increase in the maximum density. We can identify this 
phase in the evolution of the cloud as one in which turbu- 
lence is more important than self-gravity for determining the 
peak density. Although dense substructure exists within the 



cloud, created by the combination of thermal instability and 
turbulent compressions, the majority of this substructure is 
not gravitationally bound (see, for example, the discussion 



Banerjee et al.||2009| ). 



This phase in the cloud's evolution comes to an end 
once gravitationally bound cores form and begin to undergo 
runaway gravitational collapse. This occurs at t ^ 4 Myr for 
the fast flow and t 16 Myr for the slow flow, and is fol- 
lowed very quickly by the onset of star formation within the 
clouds. At the point at which the first gravitationally-bound 
core goes into runaway collapse, the minimum gas temper- 
ature has dropped to around 10 K, typical of the tempera- 
tures observed in nearby, star- forming cores (see e.g. Bergin| 



Tafalla|2007 ). Comparing the densities and temperatures 
in Figure [3] to those in Figure [5] we see that the low temper- 
atures are indeed found within the dense, collapsing cores. 
The lack of scatter in the density-temperature relation above 
10^ cm~^ reflects the fact that the thermal coupling between 
the gas and the dust becomes effective at these densities 
within prestellar cores. 

A more detailed picture of the distribution of densi- 
ties and temperatures in the flow is given in Figure |4] In 
the left-hand panels, we show how the mass of gas denser 
than a chosen density threshold nthr evolves as a function 
of time for values of nthr ranging from 10 cm~^ to 10^ cm~^ 
for the two flows. In the right-hand panels of Figure [4j we 
look at the temperature distribution of the gas. For a static 
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gas distribution in pressure equilibrium with solar metal- 
licity and with our chosen environmental parameters (ra- 
diation field strength, cosmic ray ionization rate, etc.), two 
stable thermal phases exist: a warm phase with temperature 
T > 8500 K and a cold phase with T < 350 K. Gas in the 
temperature range 350 < T < 8500 K is not thermally sta- 
ble, and will either heat up until it joins the warm phase, 
or cool down until it joins the cold phase, depending on its 
starting density. In Figure [4] we show how the gas mass in 
these three regimes evolves with time in the two flows. We 
also track the mass of very cold gas (T < 30 K). Most of 
the CO emission from local molecular clouds comes from 
gas with temperatures below around 20-30 K, and this very 
cold phase is also the most likely site of star formation in 
the cloud (see e.g. |Krumholz, Leroy & McKee 2011). 

Despite this similarity, many other features of the den- 
sity and temperature are seen to be very different when one 
compares the two flows. For example, we see that the fast 
flow maintains a much higher fraction of hot gas than the 
slow flow, reflecting the stronger shock-heating from the col- 
lision of the streams. Further, we see that the fast flow gen- 
erates much more very cold (T < 30 K) gas, and that this 
cold gas makes its flrst appearance much earlier in the simu- 
lation. By the onset of star formation, the fast flow has over 
8000 M0 of gas in the very cold state, while the slow flow has 
only around 3000 Mq in the same state. This means that 
for the fast flow, roughly 50% of all of the gas colder than 
350 K is in this very cold state, whereas the corresponding 
flgure for the slow flow is only 10%. On the other hand, the 
fast flow also has much more gas still in the unstable regime 
(350 < T < 8500 K), owing to the fact that in this case, star 
formation begins before all of the low density, unshocked gas 
has had time to cool sufliciently to reach the cold phase. 

As well as forming signiflcantly more cold gas, we also 
see that the fast flow forms more dense gas. Indeed, gas 
at number densities n ^ 10^ cm"*^ - the minimum density 
characteristic of the structures identifled observationally as 
pre- stellar cores - appears after only about 2.5 Myrs in the 
fast flow, in comparison to 19 Myr in the slow flow. We also 
see that the evolution of the density fractions in the slow 
flow is more stochastic, indicating that much of the structure 
that is formed during the collision between the streams is 
transient. This transient structure is either ripped apart by 
the ram pressure of the flows or pushed apart by internal 
thermal pressure and turbulent motions. In the fast flow, 
we do not see this behaviour; typically, the mass above each 
of our threshold densities is continually rising. 



4 CHEMICAL AND OBSERVATIONAL 
TIMESCALES 

In this section, we flrst give an overview of the general chem- 
ical evolution of the flows, and how long it takes to form 
a 'molecular' cloud in each case - that is, one that would 
be seen by an observer via CO emission. We then go on 
to look at how the post-processed CO maps of |Heitsch fc| 
Hartmann (2008) compare to our fully self-consistent and 



time-dependent treatment of the cloud chemistry. Finally, 
we look at how the observable properties of the CO vary 
with time as the clouds (and star-forming regions) are as- 
sembled. 
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Figure 5. The gas temperature-density distribution in the flows 
at the onset of star formation. 



4.1 General chemical evolution of the flows 

An overview of the chemical state of the gas can be found in 
Figure [6] The left-hand plots depict how the global chemical 
state of the gas evolves as the flow advances. They show 
the fraction of the available hydrogen that is in the form 
of H2, and the fraction of the available carbon that is in 
the form of C^, C, or CO. The fraction of the total carbon 
that is incorporated into other molecules, such as HCO^, is 
always very small and is not plotted. The right-hand plots 
show the maximum abundances of H2 and CO within the 
simulation, which tells us whether there are any molecular- 
dominated regions within the flow. Note that in this plot, 
the abundances are given with respect to the overall number 
of hydrogen nuclei (a conserved quantity), such that gas in 
the form of pure H2 will have a fractional abundance of H2 
that is 0.5, whereas gas in which all of the carbon is in the 
form of CO will have a fractional abundance of CO that is 
1.4 X 10"^. 

We start by looking at the evolution of the H2 in the 
cloud in Figure [6] The left-hand plot shows that the gas goes 
from being completely atomic - as in our initial conditions 
- to having around 10 percent of its hydrogen in molecular 
form by the point at which star formation sets in (^ 7 per- 
cent in the case of the slow flow, and ^12 percent in the 
case of the fast flow) . The initial rise in the amount of H2 is 
also sharp, going from essentially zero to around a percent 
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over a period of less than 2 Myr in each flow. Such a rapid 
rise can be understood by looking at the density evolution in 
Figure [4] We see that for each flow, the sudden rise in the H2 
fraction is accompanied by a rapid rise in the amount of gas 
with a density above 100 cm~^. Sin ce the formation time 
of H2 is of the order of 10^ /n Myr (Hollenbach McKee 



|197 9), where n is the number density of the gas, we see that 
once the gas density exceeds 100 cm~^, the time required to 
convert a large fraction of the hydrogen to molecular form 
becomes of the order of a few Myr. Therefore, the sudden 
appearance of H2 is simply a consequence of the structure 
that is formed in the flows. 

Figure |6] also shows that some pockets of gas can be- 
come almost fully molecular very early in the calculation, as 
shown the right-hand plots of Figure |6] Again this is simply 
a reflection of the density evolution that we see in Figure 
|4] What is interesting is that these pockets of H2 appear 
very early in the flows' evolution, well before the onset of 
star formation. In the slow flow, the pockets of molecular 
gas appear more than 10 Myr before the first star- forming 
core, and even in the fast flow, these pockets precede the 
star formation by about 3 Myr. In both cases, the regions of 
molecular hydrogen appear long before the flows have been 
able to assemble anything that could be construed as a star- 
forming cloud. As such, they exist in relative isolation during 
the pre-assembly phase. 

In contrast to the early appearance of H2 , the appear- 
ance of CO occurs extremely late. The simulations start with 
all of their carbon in the form of C^, and we see that most 
of the carbon stays in this form as the flows evolve, even 
once star formation has begun within the dense gas. Very 
shortly after the start of the simulation, a small fraction of 
the recombines, yielding neutral atomic carbon, and the 
amount of this that is present in the flow rises over time as 
the amount of cold, dense gas increases. In particular, the 
same increase in density that is seen in Figure [3] also causes 
the amount of neutral carbon in the flows to rise by an order 
of magnitude within only 1-2 Myrs. 

The same change in the cloud's structure that causes 
this increase in the abundance of neutral carbon and that is 
responsible for the sharp rise in the H2 fraction also prompts 
a rise in the overall CO fraction (left-hand panels of Fig- 
ure [g]). However, the CO abundance still remains extremely 
small in comparison to the C or C^ abundances: less than 
one-millionth of the total amount of carbon has been con- 
verted into CO at this point in the evolution of the flows. 
From this point on, the CO fraction steadily increases un- 
til the onset of star formation. However, even then, CO is 
about an order of magnitude less abundant than neutral 
atomic carbon. 

The evolution of the peak CO abundance in the cloud 
also differs markedly from the peak H2 fraction. We see 
(from the right-hand panels in Figure[6]), that the first region 
with a CO fraction high enough to be regarded as 'molecu- 
lar', only appears in the final Myr before the onset of star 
formation. If we look back to Figures [3] and [4] we see that in 
each of the flows, the sharp rise in the CO abundance occurs 
at the same time as the appearance of gas with a number 
density n > lO^cm"^ - the density typically associated with 
pre- stellar cores. Looking more closely at the spatial distri- 
bution of CO, we find that this is not a coincidence: the 
regions in which the carbon is almost completely converted 



into CO are indeed the same self-gravitating cores that then 
collapse to form stars (as will be discussed later in Section 

We thus see that the question of when we first form a 
'molecular' cloud depends on which molecule we are talk- 
ing about. Regions with high molecular hydrogen fractions 
are produced early in the evolution of the flows and can po- 
tentially persist for a long time prior to the onset of star 
formation if the build-up of the cloud occurs slowly. On the 
other hand, high CO fractions are produced only within self- 
gravitating gas, and hence occur only around 1 Myr before 
the onset of star formation, regardless of the details of the 
flow. Most importantly, the gas is always 'fully molecular' 
by the time it forms stars. 



4.2 Comparison to post-processed CO abundances 

Prior to this work, the only numerical study of the formation 
of CO in three-dimensional colliding flow models was that 
carried out by HH08, in which they used a post-processing 
treatment to determine which regions of their simulated flow 
would have high CO fractions. They assumed that in order 
for gas to have a high CO fraction, it would have to sat- 
isfy two criteria: it must have a temperature below 50 K 
and a mean visual extinction Ay, mean ^ 1- They found that 
in their models, CO formation was triggered primarily by 
the global gravitational collapse of the cloud in a direction 
perpendicular to the inflow, and that this resulted in rela- 
tively large regions of the cloud developing high CO abun- 
dances. They also found that these high CO abundances 
preceded the onset of star formation by around 4-5 Myr. In 
this section, we examine how the CO fractions predicted by 
their post-processing algorithm compare with those we ob- 
tain from our time-dependent chemical network, and thereby 
test whether the assumptions that HH08 use in their study 
are valid. 

In our calculations, the visual extinction is calculated 
during the tree-walk to get the gravitational forces, using 



our recently developed TreeCol algorithm (Clark, Glover 



& Klessen 2012), as described in Section 2.1 This yields 



a 48 pixel map of the column densities (and hence visual 
extinctions) seen by each SPH particle. We use this map to 
compute a mean visual extinction Ay, mean for each particle, 
using the following expression: 



2.5 



In 



48 

-y 

48 ^ 



exp(-2.5Av,i) 



(1) 



where Av,i is the visual extinction associated with pixel 
number i, and we sum over all 48 pixels. The weighted mean 
that we calculate in this fashion accounts for the fact that 
the photodissociation rates of molecules such as CO and the 
photoelectric heating rate of the gas all depend on exponen- 
tial functions of Ay, rather than directly on Ay itself. 

In Figure [7| we show the distributions of temperature 
and mean visual extinction for the particles in both flows, 
at point 0.8 Myr after the onset of star formation (recall that 
star formation occurs at a time of 16 Myr in the slow flow, 
and 4.4 Myr in the fast flow). The points are coloured by 
the CO abundance of the corresponding SPH particle, and 
results are only shown for SPH particles with fractional CO 
abundances greater than 10~^. This means that very little 
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Figure 6. Chemical evolution of the gas in the flow. In the left-hand column, we show the time evolution of the fraction of the total 
mass of hydrogen that is in the form of H2 (red solid line) for the 6.8 km s"-*^ flow (upper panel) and the 13.6 km s""*^ flow (lower panel). 
We also show the time evolution of the fraction of the total mass of carbon that is in the form of C+ (green dashed line), C (orange 
dash-dotted line) and CO (blue dot-dot-dot-dashed line). In the right-hand column, we show the peak values of the fractional abundances 
of H2 and CO. These are computed relative to the total number of hydrogen nuclei, and so the maximum fractional abundances of H2 
and CO are 0.5 and 1.4 x 10~^ respectively. Again, we show results for the 6.8 km s"-*^ flow in the upper panel and the 13.6 km s"-*^ flow 
in the lower panel. Note that the scale of the horizontal axis differs between the upper and lower panels. 



of the warm gas in the flows appears in the plot, since this 
material typically has a CO abundance below this threshold. 
In the case of the slow flow, we see from the figure that high 
CO fractions are present only in gas with a mean visual 
extinction Ay, mean > 2 and a temperature lower than 20 K. 
In this case, adopting the higher temperature threshold and 
lower extinction threshold used by HH08 would lead to an 
overestimate of the the fraction of the flow that would be 
traced by CO, and thus identified as a molecular cloud. In 
the fast flow, the HH08 approach fares better. In this case, 
the transition from low to high CO abundances does indeed 
occur for Ay, mean ^ 1, although for gas with a temperature 
below around 20-30 K; gas with T 50 K and Ay, mean ^ 1 
has a very low CO abundance. 

The motivation for the extinction threshold adopted in 
HH08 was provided by the study of |van Dishoeck Black] 
(1988), who demonstrated that in uniform density, plane- 



parallel cloud models, CO is abundant in regions that have 
visual extinctions Ay > 1. However, the [van Dishoeck &:"] 
Black (1988) study did not account for density inhomo- 



geneities in the gas, which have the effect of significantly 
complicating the relationship between the CO abundance 
and the visual extinction (see e.g. Glover et al.||2010 ). The 



HH08 approach therefore implicitly assumes that all of the 
gas with Av,mean ^ 1 and T < 50 K is dense enough to 
sustain a high CO abundance. As we see in Figure [7l this 
is not the case in our model clouds: the scatter in the rela- 
tionship between Ay, mean and n is quite large, and so much 
of the gas with Ay, mean 1 has a density below 1000 cm~^ 
and hence although it is relatively cold (T < 50 K), it is 
nevertheless too diffuse to have a high CO abundance. A 
similar density threshold of around 1000 cm~^ is also found 



in more idealised turbulent-cloud models (see e.g. [Molina 
eraL][20TT| . 



One can see by just how much the true mass of the 
clouds would be overestimated by looking at Figure [8] which 
shows the cumulative temperature distribution in terms 
of the flow mass. We see that there is around an order 
of magnitude more gas with a temperature in the range 
20 < T < 50 K than has a temperature T < 20 K. Given 
that the mean visual extinction of gas with T < 50 K is gen- 
erally quite high, as shown in Figure [71 we see that the HH08 
approach will typically result in inferred 'molecular clouds' 
that are about 10 times more massive than those predicted 
by our more self-consistent approach. 

Nevertheless, it is clear that the basic motivation under- 
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Figure 7. Weighted mean of the visual extinction seen by each SPH particle, plotted as a function of temperature (left), and density 
(right). The colour denotes the fractional abundance of CO, relative to the number of hydrogen nuclei. An SPH particle with all its 
carbon in the form of CO, would thus have an abundance of 1.4 x lO"'^. Note that most of the cloud has a CO abundance significantly 
below 1 X 10~^, and so is not visible on this plot, due to the scaling. In these plots we have taken the particle data at a point 0.8 Myr 
after the onset of the star formation, to give the cloud time to assemble. The visual extinction for each particle is computed via the 
column densities obtained by the treecol algorithm (Clark et al. 2012), as described in the text. 



lying the HH08 approach is sound. Our results suggest that 
instead of using a threshold of T < 50 K and Ay, mean > 1, 
the values T < 20 K and Ay, mean > 2 would give a more ac- 
curate picture of the location of the CO-bright gas, at least 
in the early stages of cloud assembly that we focus on here. 



4.3 Detection via CO emission 

The plots in Figure[6]show us that CO becomes highly abun- 
dant in some regions of the cloud shortly before the onset 
of star formation. However, this on its own does not tell us 
whether these regions would actually be detectable in CO 
emission. In this section, we post-process our results using 
the RADMC-SEE] radiative transfer code in order to con- 
struct maps of the ^^CO (1-0) line emission coming from the 
region immediately surrounding the first star-forming core 
in each flow. To compute the CO level populations, we use 



the large velocity gradient approximation (Sobolev 1957), 
as implemented in RADMC-3D by |Shetty et al.| ( |2011a[b ). 
The SPH data is interpolated onto a 256^ grid of size 6.48 
pc, centred on the star-forming core. In this manner, the 



^ http:/ /www. ita.uni-heidelberg.de/~dullemond/software/radmc- 
3d/ 



self-consistently computed abundances of the CO from the 
simulations can be used as the basis of the radiative transfer 
calculation. The radiative transfer is done along the x-axis 
of the grid, such that the rays are fired from negative x to 
positive X. The resulting maps are hence those that would 
be seen by an observer sitting at positive x (i.e. one who is 
looking along the flow). 

The images in Figure [9] show the results from the ra- 
diative transfer calculations, along with the column number 
density for the same region. The line transfer data is shown 
as the velocity-integrated intensity {Wco) - the quantity 
that is commonly used to determine H2 column densities 
via the so-called 'X-factor' ( Solomon et al.p987 Young & 
'Scoville| |1991 Dame et al.||2001 ). We perform the radiative 
transfer at 4 times in each flow's history, spanning a period 
from 2 Myr before the onset of star formation, to 0.8 Myr af- 
ter (which is also the point at which we halt the simulations 
due to the rapidly increasing computational expense). 

We see that in both cases, Wco is large in the region 
immediately surrounding the star-forming core, and that the 
dense prestellar core is embedded in a larger region of diffuse 
CO emission. The range of values for Wco in this region - 
from a few K kms~^ to around 20-30 K kms~^ - appears to 
be consistent with the range of values found in local clouds 



© 0000 RAS, MNRAS 000, 000-000 



Molecular cloud formation timescales 11 



10^ 



100 



10 



1 = 




1 10 100 10-" 

temperature [K] 

Figure 8. The cumulative temperature distribution in the slow 
and fast flows at the onset of star formation. 



such as the Taurus molecular cloud ( Goldsmith et al.||2008 



Narayanan et al. 2008). However, the length-scale associ- 
ated with the CO-bright regions does differ between the two 
flows. For the slow flow, the region around the prestellar 
core is 2 pc in length at the onset of star formation, while in 
the fast flow the region is somewhat smaller, having a diam- 
eter of roughly 0.5 pc. In both cases, the extent of the CO 
emission reflects the column density distribution in the gas. 
In general, however, it seems that at these early times the 
CO is confined to regions that are undergoing gravitational 
collapse. 

The images shown in Figure [9] also demonstrate that 
both the size of the CO-bright region and the strength of 
the emission from this region vary strongly with time. At 
a time of 2 Myr before the onset of star formation, there 
is essentially no visible CO emission coming from the fast 
flow, and only one small region of relatively faint emission 
in the slow flow|^ However, the size of the CO-bright re- 
gion and the strength of the emission from this region both 
increase rapidly as the prestellar core begins to undergo 
gravitational collapse. From the column density images, we 
can see that this change largely just reflects the change in 
the column density distribution that occurs as the core col- 
lapses. Although the observational definition of 'detection' 
of a molecular cloud via CO emission can vary, depending 
on the method used to extract the cloud from the data- 
set and the telescope used to conduct the survey, a value of 
Wco ^ IKkms"^ appears to be a practical threshold for the 
robust detection of 'molecular gas' (see, for example, the dis- 
Goldsmith et al.|2008| and |Heyer et al.|2009| ). In 



cussions m 



our maps, we see that Wco falls off rapidly below this value, 
and so we can conclude that the emission seen in Figure [9] 
is detectable. The Wco maps therefore confirm the picture 
that we discussed in Section [4. 1| namely that the emergence 
of CO occurs fairly late in the evolution of the flow, and at 
roughly 2 Myr before the onset of star formation. 



^ Strictly speaking, there is a very low level of CO emission com- 
ing from every part of the cloud, but in the image we only indicate 
the emission when V^co ^ 0.1 K km s"-*^, as emission below this 
threshold would not in practice be detectable. 



From this we can conclude that the colliding flow model 
can indeed produce an observationally detectable molecular 
cloud prior to the point at which star formation begins in 
the cloud. Furthermore, the timescale over which the CO 
becomes observable is typically around 2 Myr, and appears 
to have little dependence on the speed of the flow that is 
forming the cloud. 



5 THE ONSET OF STAR FORMATION 

Our discussion so far has touched upon conditions in the 
clouds that accompany the onset of star formation. In this 
section we take a closer look at how the gravitational insta- 
bility occurs and compare our results with previous studies 
of self- gravitating colliding flows. 

The plots in Figure [To] show how the properties of the 
gas vary as a function of distance from the first collapsing 
core to form in each simulation. We plot the radially aver- 
aged number density and temperature of the gas, along with 
the mass enclosed within a sphere of that radius. We also 
examine how the ratios Ec/Et and Eg/{Ek + ^t) vary 
with distance from the core, where Eg is the magnitude of 
the gravitational energy, Et is the thermal energy and Ek 
is the kinetic energy. Results for the slow flow are shown on 
the left, and those for the fast flow are shown on the right. 
The profiles are taken from the last output snapshot pro- 
duced by each simulation before the onset of star formation, 
and hence correspond to a time less than 10^ years prior to 
the formation of the first sink particle]^ 

As in the previous sections, we again see a very differ- 
ent picture when we compare the conditions surrounding the 
star-forming core in each flow. The main difference between 
the simulations is that the slow flow is almost entirely Jeans 
unstable by the onset of star formation (that is, its grav- 
itational energy is greater than its thermal energy), while 
the fast flow is in general gravitationally stable. This is due 
primarily to the fact that turbulence in the fast flow has not 
had time to trigger the cooling instability, and so most of 
the gas is still in the warm phase. This is apparent when 
we compare the radial temperature profiles in the each flow; 
the temperature rises gradually in the slow flow, while in the 
fast flow we see a sudden rise in the temperature at around 
1 pc from the centre of the core. This sudden temperature 
rise marks the boundary of the collapsing region in the fast 
flow. In the slow flow, we see instead that this boundary 
is marked by a rise in the kinetic energy, indicating that it 
is the turbulent motions in the flow that determine which 
regions can locally collapse and which cannot. 

The fact that the entire cloud becomes gravitationally 
unstable in the slow flow is broadly consistent with the 
previous numerical studies of self-gravitating colliding flows 



(|Vazquez-Semadeni et al.|2006||Heitsch et al. 


2006| Vazquez- 


Semadeni et al.|2007| 


Hennebelle et al.|2008 


|Banerjee et al. 


2009 


[Vazquez- Semac 


eni et al.||2009| Rosas-Guevara et al. 


2010 


). Typically it is found that the first star- forming cores 



appear while the post-shock layer itself is undergoing col- 
lapse (see the discussions in|Heitsch, Hartmann & Burkert 



2008 and Vazquez- Semadeni et al. 2009). In contrast, the 



We produce output snapshots every 10^ yr during the runs. 
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behaviour seen in our fast flow seems to differ from the stan- 
dard picture, in that the collapsing regions are localised, and 
not accompanied by the collapse of the post-shock layer. The 
origin of this difference is unclear, and the fact that the ex- 
isting studies (at least the high resolution, self- gravitating 
studies) do not present data similar to that expressed in our 
Figure |4] means that a detailed investigation of this differ- 
ence lies beyond the scope of this paper. However one poten- 
tial difference between our results and the previous studies 
is that we include a more self-consistent treatment of the 
photoelectric emission heating, as we are able to accurately 
calculate the attenuation of the external radiation field via 
our TreeCol algorithm. Thus it is likely that our fraction 
of very cold gas (i.e. gas colder than 30 K) is much higher, al- 
lowing gravitationally bound structures to form more readily 
in our simulations. 

A final interesting feature regarding the onset of star 
formation in our simulations is the fact that both flows have 
roughly the same mass in the local collapsing region - at 
least to within a factor of two. Given the very different con- 
ditions in the flow at the point where the gravitational insta- 
bility sets in, this is perhaps somewhat surprising. However, 
the temperature of the gas in the cold regime is a strong 
function of density, as can been seen in Figure [5] Although 
the fast flow does have slightly higher temperatures on aver- 
age at a number density of 10^ cm~^ - the density at which 
the gravitational instability sets in, as shown in Figure [lO] 
- the difference is fairly small. As such, the Jeans mass in 
the collapsing regions is roughly the same. This ability of 
the gas to regulate its temperature is obviously crucial for 
explaining why the conditions in regions of nearby star for- 
mation are so similar ( Larson||l985| 2005 Elmegreen et al 
20081. 



and forms a region that grows until it is gravitationally un- 
stable. This behaviour is similar to that seen in the simula- 



6 DISCUSSION 

The results of this study show that it is possible to assem- 
ble star-forming regions in two ways - either by putting to- 
gether pre-existing, cold molecular clumps, or forming them 
directly from the warm ISM via the cooling instability. The 
former scenario is what we see in our 'slow' calculation, in 
which the majority of the gas in the incoming flow has time 
to cool before hitting the shock at the edge of the central 
cloud. As such, the gas is cold and clumpy by the time it en- 
counters material coming from the other direction, and thus 
the set-up is more akin to the well-studied 'cloud-cloud' col- 
lision scenario rather than the classic colliding flow set-up, 
in which only the interface between the flows has been per- 
turbed. In fact, our set-up is similar to that employed by 



Rosas-Guevara et al. (2010), in which they show that the 
exact properties of the star formation in the clouds are sen- 
sitive to the initial level of turbulence that is injected into 
the colliding flows. 

The second scenario, in which the dense bound regions 
form directly from the cooling instability in the shocked re- 
gion, is represented by our 'fast' flow calculation. In this 
case, although the flow is perturbed throughout, there is 
insufficient time for the thermal instability to occur in the 
low density inflowing gas and so the majority of the gas 
entering the shock is still in the unstable or warm phases 
(T > 1000 K). In this case, the warm, incoming gas cools 



tions of Banerjee et al. (2009), although as discussed above 
in Section |5] we do not see the accompanying collapse of 
the post-shock layer in this case. As the temperature in the 
post-shock gas is set by the balance between photo-electric 
emission and line-cooling, it is similar to that found in 
the star- forming region that forms in the slow flow. There- 
fore, although the assembly process is very different in the 
two cases, the resulting star-forming regions have similar 
properties. 

Which of these two scenarios is more likely to occur 
in reality? [Heitsch et al. ( 2006 ) point out that the density 
regime from 1 to 10 cm~ is thermally unstable, and so it is 
likely that clouds are at least in part assembled via collisions 
between cold gas. The conditions that we see in our fast flow 
are therefore likely to be somewhat artificial, as in reality the 
flows would have had ample time to undergo the thermal 
instability. Our calculations also suggest that much of the 
cold gas that participates in the cloud formation process 
may already consist of H2 long before the final 'cloud' (or 
star-forming region) is assembled. This gas would not be 
observable via its CO emission until it has become part of 
the main cloud. As such, our simulations support the picture 
outlined by Pr ingle, Allen Lubow] ( [200 1), who propose 
that clouds form out of pre-existing, CO-dark, molecular 
gas. Similar behaviour was also suggested by the study by 
[Dobbs et al.| ( |2008| ), which looked at GMC formation in 
spiral arms. In this context, it is interesting to note that 
there is growing observational evidence for the existence of 
significant quantities of H2-rich, CO-poor gas in the Milky 



Way ( 


Grenier et al.||2005 


Langer et al.||2010 Pineda et al. 


2010 Velusamy et al. 2010 Ade et al.|2011 ). Our simulations 



suggest that this so-called "dark" molecular gas is a natural 
outcome of the process of cloud assembly. 

Finally, we note that as the CO appears fairly late in 
the cloud assembly process, it plays only a minor role in 
the thermodynamics of cloud formation. Throughout much 
of the assembly of the cloud, the cooling is dominated by 
atomic line cooling from C^ and O, with CO becoming im- 
portant only in regions where the gas is already cold and 
dense. This is consistent with the picture presented in | Glover] 
Clark|p012 a), where it was argued that molecules are not 



needed to form stars, but rather accompany the star forma- 
tion process. 



7 CONCLUSIONS 

We have presented results from the first three-dimensional 
calculations to follow the non-equilibrium chemistry of H2 
and CO formation from the warm neutral medium to the 
cold molecular regime, in the context of the colliding flow 
picture of molecular cloud formation ( Elmegreen|199 3^ . Our 
time-dependent chemical network is coupled to a detailed 
model of the thermodynamics of the ISM that accounts for 
the main heating and cooling processes, such as atomic and 
molecular line cooling (from H, C^, C, O, CO and H2), 
photoelectric emission and cosmic ray heating, and the en- 
ergy transfer between the gas and dust. The chemical and 
thermodynamical model are incorporated into the publicly 
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Figure 9. The images show the evolution of the column number density, N , and the velocity-integrated intensity in the J = 1 — ^ line 
of ■'^■^CO, Wco (1-0), for the region in which the first star forms in each of the fiows. Four times are shown: 2 Myr prior to star formation 
(upper left-hand panels), 1 Myr prior to star formation (upper right-hand panels), the point of star formation (lower left-hand panels), 
and 0.8 Myr after the onset of star formation (lower right-hand panels). The CO integrated intensity map is obtained via a radiative 
transfer calculation performed with the RADMC-3D code, and uses the large velocity gradient approximation to compute the CO level 
populations. 



available smoothed particle hydrodynamics code Gadget2 
( Springel|20"05t . 

We have presented calculations in this paper that fol- 
low the evolution of two colliding flows with velocities 6.8 
km s~^ and 13.6 km s~^, which we refer to as our 'slow' 



and 'fast' flows, respectively. The flows start essentially fully 
atomic in composition, except for a small ionisation frac- 
tion of 6.5 X 10~^, set by the equilibrium between cosmic 
ray ionisation and radiative recombination. The density and 
temperature in both flows is initially 1 cm""^ and 5000 K, 
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Figure 10. The panels show the radial profiles centered on the first SPH particle to find itself in a star in the 6.8 km s~^ fiow simulation 
(left) and in the 13.6 km s"-*^ fiow (right). The profiles are computed just before the onset of star formation in each case. Temperature 
and density are mass-weighted averages, while 'mass' denotes the mass enclosed at a given radius. In the bottom panel, Eq^ Ej^ and E^ 
denote the magnitude of the gravitational energy, the kinetic energy, and the thermal energy, respectively. 



respectively. To seed the structure that will eventually de- 
velop into molecular clouds, we impose a turbulent field 
with Vrms = 0.2 Cs, whcrc the isothermal sound speed is 5.6 



km 



The velocities follow the observed scaling law of 



V oc L°-^ (e.g. |Heyer k Brunt||2004 ) 



Star formation - as captured by the formation of 'sink 
particles' ( Bate, Bonnell Price|1995 ) - occurs in the calcu- 
lations at around 16 Myr for the slow flow, and 4.4 Myr for 
the fast flow. By this point, each flow has over 1000 Mq of 
gas sitting below a temperature of 30 K. We find that in 
both cases, the flows are able to form pockets of gas that 
have very high H2 fractions, and that this occurs early in 
the evolution of the flow, at a time of around 6.5 Myr in the 
slow flow, and 1.5 Myr in the fast flow. 

However the formation of regions with a high CO abun- 
dance is significantly delayed with respect to the molecular 
hydrogen. In both flows, we find that significant quantities 
of CO do not form until around 2 Myr before the onset of 
star formation. Using the results from our SPH simulation 
as input to line radiative transfer calculations, we were able 
to produce synthetic emission maps for the CO (1 — 0) line. 
These synthetic emission maps confirm that the appearance 
of CO-bright regions in the clouds precedes the formation of 
stars by no more than about 2 Myr. 

Our results suggest that even in cases when the assem- 
bly of the cloud takes place over many millions of years, 



the interval between the cloud becoming observationally de- 
tectable in CO and the onset of star formation is short, 
of the order of 1-2 Myr. This is consistent with the idea 
that molecular clouds can form out of "dark" molecular gas 
(|Wolfire et al.||2010| - cold gas that has a high H2 content 



but contains little or no CO (Pringle, Allen k Lubow|2001 



Dobbs et al.||2008| |Dobbs||2008 ) 
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